Identification of three elevenin receptors and roles of elevenin disulfide bond and residues in receptor activation in Aplysia californica

Neuropeptides are ubiquitous intercellular signaling molecules in the CNS and play diverse roles in modulating physiological functions by acting on specific G-protein coupled receptors (GPCRs). Among them, the elevenin signaling system is now believed to be present primarily in protostomes. Although elevenin was first identified from the L11 neuron of the abdominal ganglion in mollusc Aplysia californica, no receptors have been described in Aplysia, nor in any other molluscs. Here, using two elevenin receptors in annelid Platynereis dumerilii, we found three putative elevenin GPCRs in Aplysia. We cloned the three receptors and tentatively named them apElevR1, apElevR2, and apElevR3. Using an inositol monophosphate (IP1) accumulation assay, we demonstrated that Aplysia elevenin with the disulfide bond activated the three putative receptors with low EC50 values (ranging from 1.2 to 25 nM), supporting that they are true receptors for elevenin. In contrast, elevenin without the disulfide bond could not activate the receptors, indicating that the disulfide bond is required for receptor activity. Using alanine substitution of individual conserved residues other than the two cysteines, we showed that these residues appear to be critical to receptor activity, and the three different receptors had different sensitivities to the single residue substitution. Finally, we examined the roles of those residues outside the disulfide bond ring by removing these residues and found that they also appeared to be important to receptor activity. Thus, our study provides an important basis for further study of the functions of elevenin and its receptors in Aplysia and other molluscs.


Results
Molecular characterization and phylogenetic analysis of Elevenin. In a previous study 11 , an mRNA with a length of 1112 bp (NCBI accession number: k02458.1) has been cloned from L11 neurons in Aplysia, and it contains a complete CDS sequence for a precursor for a 151 aa neuropeptide (accession number: AAA27762.1) ( Fig. 2A). A mRNA sequence (TRINITY_DN18985_c0_g1_i2) was also found in AplysiaTools database (Fig. 2B), which differs from k02458.1 in one nucleotide, but the protein it encodes is identical to AAA27762.1.
We used NeuroPred 63 to predict the possible elevenin peptide from the Aplysia elevenin precursor (Fig. 2C). One putative mature peptide with 20 amino acids was produced: RPRIDCTRFVFAPACRGVSA, which contains two cysteines separated by 8 amino acids. We compared the Aplysia elevenin peptide with elevenins in other molluscs, annelids, nematodes, and arthropods ( Fig. 1, see Supplementary Table 1 for information on these sequences). Cys, Gly, Arg, Pro, and Asp residues appear to be highly conserved (Fig. 1A), whereas the total number of residues in elevenins between different species varied from 19 to 27. In addition, we generated a frequency plot (Fig. 1B) of single residues in the peptides (http:// weblo go. berke ley. edu/ logo. cgi). Both the alignment results (Fig. 1A) and the frequency plot (Fig. 1B) showed the same conserved residues. We verified the expression of the mRNA encoding elevenin precursor by performing RT-PCR on the abdominal ganglion from Aplysia (Fig. 3A, see Supplementary Fig. 1A for the complete gel, primer sequences are in Supplementary Table 2). As expected, a 456 bp fragment (identical to TRINITY_DN18985_c0_g1_i2 in the AplysiaTools database) encoding the 151 aa precursor could be amplified. We also tried other ganglions, including cerebral, buccal, and pedal ganglions, but could not obtain this sequence, supporting that elevenin mRNA may be primarily expressed in the abdominal ganglion.
Next, we compared identities and similarities between the Aplysia elevenin precursor and other elevenin precursors in selected invertebrate species (Table 1, see Supplementary Table 3 for information on these sequences). The Aplysia elevenin precursor is most closely related to that of Platynereis with a similarity of 41.52%. Furthermore, we compared 7 elevenin precursors in molluscs with 3 in annelids with BioEdit (Fig. 4), which showed that the regions for the elevenin peptide sequences are the only portions that are highly conserved and indicated that no other peptide is encoded by these precursors.

Identification of three putative elevenin receptors in Aplysia.
Currently, there are no complete elevenin receptors that have been identified in any mollusc. To identify putative elevenin receptors in Aplysia, we used the sequences of the two deorphanized elevenin receptors from the annelid Platyneris dumerilii (the first identified elevenin receptors) to perform a BLAST search of NCBI GenBank. A total of 4 Aplysia sequences (XM_005096018.3, XM_035970075.1, XM_013090833.1, XM_035970908.1) were found. The sequence (XM_013090833.1) is an incomplete sequence for the Aplysia leucokinin-like receptor (ALKR) that  61 . The three Aplysia NCBI mRNA sequences were incomplete, but we found complete sequences for all three in the AplysiaTools database (Supplemental Table 4). The pre-  www.nature.com/scientificreports/ dicted intron-exon structures of two of the receptors were shown in Fig. 5, whereas no genomic information is available for the third receptor in the AplysiaTools database because the database is still evolving and might not be complete at this point. The 7TM and conserved motifs of the three receptors were analyzed using the NCBI Conserved Domain Database (https:// www. ncbi. nlm. nih. gov/ Struc ture/ cdd/ wrpsb. cgi) and TMHMM 2.0 (https:// dtu. biolib. com/ DeepT MHMM) (Fig. 6A). The similarities between apElevR1, apElevR2, apElevR3, and the two receptors in Platynereis are approximately 54%, 55%, and 40% respectively (see the Clustal table in Supplementary Table 4), therefore we putatively named them elevenin receptors (apElevR1-3). We cloned apElevR1, apElevR2, and apElevR3 from the Aplysia CNS cDNA (  Table 2). The PCR products were found to produce identical proteins to those in the AplysiaTools. We do note that although apElevR1 and apElevR3 mRNA sequences are identical to those in AplysiaTools, apElevR2 does differ from the AplysiaTools sequence in one nucleotide, but this single nucleotide polymorphism (SNP) did not affect the protein sequence. The gene sequences have been deposited into the NCBI database (apElevR1 GenBank accession number: OQ686759; apElevR2 GenBank accession number: OQ686760; apElevR3 GenBank accession number: OQ686761). We also compared Aplysia putative elevenin receptors with elevenin receptors from several species including those from Platynereis and Ixodes using BioEdit ( Supplementary Fig. 2).
To understand the phylogenetic relationship of the three putative elevenin receptors, we first took advantage of a number of different Class-A GPCR sequences primarily from mollusc Crassostrea gigas and Lottia gigantea for a phylogenetic tree generated in our previous work (For the sequence information, see their Supplementary  Table 4) 8 . Then we added the three putative Aplysia elevenin receptor sequences, together with predicted elevenin receptors in L. giagantea and C. gigas (see Supplementary Table 4). The updated phylogenetic tree (Fig. 6B) shows that the three putative elevenin receptors of Aplysia are clustered with elevenin receptors of L. giagantea and C. gigas, supporting that they are likely elevenin receptors. Finally, we generated a phylogenetic tree with only the elevenin receptors from Aplysia, and other molluscs, annelids, and arthropods (Fig. 7). The three putative elevenin receptors in Aplysia are clustered with molluscan receptors. Activation of putative receptors by the elevenin peptide and roles of specific residues in elevenin. To determine whether the Aplysia elevenin can activate the putative elevenin receptors, we cloned apEl- www.nature.com/scientificreports/ evRs to pcDNA3.1 (+) plasmid and transfected the three apElevRs in CHO-K1 cells. We used an IP1 accumulation assay to measure the concentration of IP1, which is a by-product in the Gαq signaling pathway, to provide evidence for the activation of receptors by elevenin. Similar to previous work in Platynereis 2 , we did not need to co-transfect a promiscuous Gαq protein (also known as Gα16) to obtain an IP1 response. This indicates that the three Aplysia elevenin receptors could couple to the Gαq protein endogenously present in CHO cells. Initially, we used the elevenin peptide at two concentrations of 10 -11 M and 10 -4 M to screen apElevR1, apElevR2, and apElevR3. The results showed that 10 -4 M elevenin showed significantly higher activation on all three receptors compared with 10 -11 M, supporting that these three receptors are true receptors for elevenin (Fig. 8A). We also used synthetic elevenin peptide without the disulfide bond and found that elevenin without the disulfide bond could not activate the receptors, suggesting the importance of the disulfide bond for the activity of receptor (Fig. 8A, and E). Further, we used different concentration gradients of elevenin, ranging from 10 -12 to 10 -4 M, to generate the dose-response curves ( Fig. 8B-E). The EC 50 values of the three elevenin receptors were 25, 1.2, and 3.1 nM respectively. We have shown that C6, C15, G17, R16, P13, D5, F11, V18, and I4 in elevenins appear to be highly conserved (Fig. 1). To determine the potential importance of these conserved residues in receptor activity, we synthesized analogs by replacing these amino acids, G17, R16, P13, D5, F11, V18, and I4 (except Cys) with Ala, respectively. The results of the IP1 accumulation assay showed that the log[EC 50 ] and EC 50 values of the elevenin receptors by most of the analogs increased significantly, and the degree of effects was approximately G17 > R16 > P13 > D 5 > F11 > V18 > I4, particularly for apElevR1 (Fig. 9). Specifically, apElevR1 and apElevR3 had similar results for the peptide analogs with different conserved amino acids replaced with Ala. For apElevR1, the log[EC 50 (Table 2). In Fig. 8A, we showed that three receptors can't be activated by elevenin without a disulfide bond. We sought to investigate whether sequences other than the disulfide bond might play some roles in receptor activity. Thus, we synthesized two truncated elevenin analogs with different lengths and kept sequences between cysteines the same: Elevenin 5-17 and Elevenin 6-15 . Elevenin 5-17 has the conserved D5, R16, and G17 outside of two cysteines, whereas elevenin [6][7][8][9][10][11][12][13][14][15] only has the sequence between the cysteines. The three receptors responded differently to the two truncated elevenin peptides. The log[EC 50 ] values of the truncated elevenin 5-17 and elevenin [6][7][8][9][10][11][12][13][14][15] were all significantly increased when compared with that of elevenin for the three receptors. For apElevR1 and apElevR3, the log[EC 50 ] values of elevenin [6][7][8][9][10][11][12][13][14][15] were increased significantly compared with those for elevenin [5][6][7][8][9][10][11][12][13][14][15][16][17] (Fig. 10A,C,D,F, and G). In fact, elevenin 6-15 apparently did not have any effect on apElevR1. For apElevR2, both analogs had higher log[EC 50 ] values compared with that for elevenin ( Fig. 10B,E, and G). This is consistent with the effects of residue substitution experiments on apElevR2 (Fig. 9C,D,H and J), indicating that apElevR2 is highly sensitive to any alternations in elevenin. Thus, the amino acids outside of the disulfide bond also appear to be important for the elevenin peptide to retain its activity on the receptors.
Putative expression levels of the elevenin precursor and receptors. Finally, we characterized the expression of elevenin and its receptors using an NCBI database from a broad spectrum of RNA-seq data (GSE79231) obtained from Aplysia adult and developmental stages 64 . The highest level of the elevenin precursor (NM_001204550.1) expression was found in the CNS (Supplementary Fig. 3A). ApElevR1 (XM_005096018.2) was primarily expressed in statocyst as well as mantle and gills ( Supplementary Fig. 3B) with relatively weak expression in the CNS. ApElevR2 (XM_013083173.1) was also expressed most highly in statocyst, and with the Table 1. Identities and similarities of the elevenin precursor from Aplysia with elevenin precursors from other molluscs, annelids, nematodes, and arthropods using BioEdit (pairwise alignment-calculate identity/similarity for two sequences). www.nature.com/scientificreports/ second high-level expression in the CNS (Supplementary Fig. 3C). Note that XM_013083173.1 (apElevR2) is incomplete and has been replaced by the accession number of XM_035970075.1 in NCBI, which is still incomplete (see the section "Identification of three putative elevenin receptors in Aplysia"). In addition, no information is available on apElevR3 (XM_035970908.1).

Discussion
In previous work, elevenin was initially identified as a cDNA sequence encoding L11 neuropeptide precursor from Aplysia 11 . Here, we verified the elevenin precursor by nucleic acid sequencing, and more importantly, we identified three elevenin receptors in Aplysia for the first time and demonstrated the roles of the disulfide bond and conserved amino acid residues of the elevenin ligand on receptor activity through bioinformatics and cellbased assay.  www.nature.com/scientificreports/ Elevenin peptides and their precursors. After the elevenin precursor was first discovered in L11 neurons of Aplysia 11 , homologous precursors have subsequently been found in other molluscs and C.elegans, and therefore they were renamed elevenin 12 . Later, an elevenin precursor was identified by sequencing the cDNA in Ixodes 9 . Subsequent bioinformatic work has predicted the presence of the elevenin precursor and peptides in more species, including molluscs, annelids, nematodes, and arthropods ( Fig. 1). Consistent with the previous work, we were able to clone the Aplysia elevenin precursor from the abdominal cDNA (Fig. 3A). It encodes a 20 amino acid peptide if cleavage is done after the first Arg residue following the signal peptide. A few studies have identified the mature forms of elevenins in two species (Carabus violaceus and Cataglyphis nodus). Specifically, the elevenin peptide likely with the disulfide bond was identified by mass spectrometry (MS) in the brains of Cataglyphis nodus 18 . In addition, it was noted that the elevenin in Carabus violaceus has been identified by MS, but a clear MS spectrum was not obtained, although there were ion signals in mass fingerprints 65 . Interestingly, conopeptides in the venoms of cone snails contain peptides often with two or more cysteines 66 . Indeed, venoms of cone snails include precursors that encode peptides similar to elevenins, and the mature elevenin-like peptide with a disulfide bond in the Australian cone snail Conus victoriae, named elevenin-Vc1, was identified by MS 67 . Thus, although it remains to be determined experimentally in more species, we expect that the mature forms of endogenous elevenins in protostomes are likely to form a disulfide bond. This is supported by our finding that elevenin without the disulfide bond cannot activate the elevenin receptors at all in Aplysia (Fig. 8).
Identification of three Aplysia elevenin receptors. Although the function of elevenin was first studied in C.elegans, its actual elevenin receptors have not yet been found in this species. In fact, prior to our work, there are a total of 6 elevenin receptors that have been described (Fig. 7). First, there are two GPCRs that have been deorphanized as elevenin receptors in both annelid Platynereis and arthropod Ixodes 2,9. There is a single elevenin receptor in Nilaparvata lugens or Bombus terrestris 23  In molluscs, before this present work, there is only one study that identified a partial sequence of a putative elevenin receptor in sea slug Plakobranchus ocellatus 26 , which appeared to contain only 2 transmembrane domains based on our analysis. In our study, we identified three complete elevenin receptors in Aplysia. The EC 50 values of the three receptors apElevR1, apElevR2, and apElevR3 in response to elevenin with the disulfide bond were 25 nM, 1.2 nM, and 3.1 nM, respectively (Fig. 8), which are in the similar ranges to the EC 50 values for receptors www.nature.com/scientificreports/ www.nature.com/scientificreports/ in the Platynereis and Ixodes. Moreover, for the assessment of receptor activity of all elevenin receptors in the three species, an expression system, i.e., CHO cells, is used and a promiscuous Gαq protein (also known as Gα16) was not needed to obtain a response for IP1 (Aplysia) or Ca 2+ (Platynereis and Ixodes), both of which are generated in the Gq signaling pathway. Thus, all elevenin receptors identified so far could couple to the Gαq protein endogenously present in CHO cells. However, it remains to be determined whether these receptors could also couple to other Gα proteins, such as Gαs or Gαi. Overall, our work is the first time to explicitly identify the complete elevenin receptors in a mollusc. Given that currently, only Aplysia has three elevenin receptors, whereas other species had at most two 2,9 , an open question is whether the three elevenin receptors are only present in molluscs, or whether there may be additional receptors that remained to be discovered in other species, e.g., in annelid Platynereis. This issue could be resolved as more elevenin receptors are being characterized in protostomes in the future. It is also interesting to note that phylogenetic trees (Figs. 6 and 7) suggest that each of apElevR1, apElevR2, and apElevR3 might constitute a separate clade with putative orthologs of other molluscan species (C. gigas and L. gigantea), which differ from the two receptors in I. scapularis. The results imply that at least three elevenin receptors might have been generated in a common ancestor of molluscs, and are conserved in each current mollusc, although some of them might have been deleted or become pseudogenes during their evolutionary processes. In contrast, the two Ixodes elevenin receptors could be generated by duplication in the evolutionary lineage of Ixodes. As elevenin receptors are being identified in more species, we will have a better idea of whether these evolutionary insights hold.

Role of elevenin disulfide bond and conserved residues in receptor activation.
We showed that elevenin analogs without the disulfide bond could not activate any of the three receptors (Fig. 8), supporting that the mature elevenin in Aplysia likely forms a disulfide bond and the disulfide bond is important for receptor activity. We recently identified another Aplysia peptide with a disulfide bond, i.e., allatostatin-C (SHYSSMCM-FNVVACY) 8 . In this case, Aplysia AstC without a disulfide bond could activate its receptor: AstC-R, although with a significantly higher EC 50 value compared to that in response to AstC with the disulfide bond. It is notable that there are 6 residues between the two cysteines in AstC, compared to 8 in elevenin. Perhaps, the smaller www.nature.com/scientificreports/ number of residues between the two cysteines allowed some of the AstC to form a disulfide bond spontaneously so that AstC without the disulfide bond could activate AstC-R at a lower potency. For comparison, for the Aplysia vasotocin (a homolog of mammalian vasopressin/oxytocin) with a disulfide bond (CFIRNCPKG-amide, 4 residues between the two cysteines) we characterized recently 62 , the analog without the disulfide bond actually had similar EC 50 values on the two vasotocin receptors to the vasotocin with the disulfide bond. However, when the cysteines are protected by acetamidomethy to prevent the possible spontaneous formation of the disulfide bond, the receptor activity was dramatically reduced, supporting the importance of the vasotocin disulfide bond for receptor activity. Our alanine substitution experiments demonstrated the specific roles of conserved residues in elevenin for receptor activation. We selected seven residues that are the most conserved ones except the two cysteines, i.e., G17, R16, P13, D5, F11, V18, and I4 with their degree of conserveness in this order (Fig. 1). Indeed, for apEl-evR1 and apElevR3, the substitution of G17 with alanine had the largest effect, and substitution of I4 and V18 had lowest effects (Fig. 9, Table 2), as also demonstrated for other neuropeptides 61 . Interestingly, for apElevR2, although substitution of G17 with alanine tended to have the largest effect, the substitution of the other five residues also had similarly profound effects, indicating apElevR2 was more sensitive to residue substitution compared to apElevR1 and apElevR3. For comparison, in several other peptides with the disulfide bond, alanine substitution of residues within the disulfide bond tends to have a larger impact on receptor activity than those outside the disulfide bond 62,68,69 . Thus, different receptors might have different sensitivities to alanine substitution in peptides with a disulfide bond. www.nature.com/scientificreports/ www.nature.com/scientificreports/ Given that the most conserved residues are located between D5 and G17, we also made a truncated elevenin analog (elevenin [5][6][7][8][9][10][11][12][13][14][15][16][17] ) that included only these 13 residues, and another one (elevenin [6][7][8][9][10][11][12][13][14][15] ) that only included the two cysteines and the 8 residues between them (a total of 10 residues). If the elevenin activity depended primarily on the disulfide bond and the conserved residues, we would expect a small impact of elevenin 5-17 on receptor Table 2. The EC 50 values on apElevRs by elevenin analogs with alanine substitutions. Analogs are ordered based on the conserveness of each residue (Fig. 1). Analogs with substituted residues inside the disulfide bond were highlighted in italic.  www.nature.com/scientificreports/ activity. However, this is not what we observed (Fig. 10), indicating that the other residues and/or the overall peptide structure maintained by the complete 20 residues also played some roles in receptor activation. We did find that the shorter elevenin 6-15 had a bigger impact than elevenin [5][6][7][8][9][10][11][12][13][14][15][16][17] , which is expected because there are fewer conserved residues in elevenin [6][7][8][9][10][11][12][13][14][15] (Fig. 10). Again, apElevR2 was more sensitive to residue deletion compared to apElevR1 and apElevR3. Overall, both alanine substitution and residue deletion experiments indicated that even though the disulfide bond is essential, conserved residues, other less conserved residues, and/or the overall structure of the peptide might also play roles in receptor activation. Moreover, apElevR2 is more sensitive to alternation of the elevenin ligand than apElevR1 and apElevR3. Along this line, it is worth noting that apElevR2 also differed from apElevR1 and apElevR3 in that it had relatively low maximum IP1 responses (Fig. 8). Future work is needed to determine how the elevenin and its analogs interact with these receptors, perhaps using the approach involving AI protein structure prediction-based modeling described recently 61 . Such studies may inform us how and why the three receptors might have different sensitivities to alanine substitutions and residue deletion of the elevenin.
In summary, we have identified three elevenin receptors in Aplysia, which, together with the elevenin precursor, constitute a complete characterization of the elevenin signaling system in a mollusc. Given that several neural circuits that mediate a variety of behaviors, such as feeding and locomotion 34,36,40,70 , have been studied in Aplysia, it would be of interest to determine whether elevenin might play a modulatory role in these circuits. Indeed, RNA-seq data of different tissues (Supplementary Fig. 3) showed that the elevenin precursor is most highly expressed in the CNS, whereas apElevR1 and apElevR2 are not expressed most highly in the CNS, but they are present in the CNS. These data suggest that the elevenin signaling system might play some roles in these and/ or other neural circuits. In addition, apElevR1 and apElevR2 are also highly expressed in the statocyst, and/or the mantle and gills, suggesting that the elevenin signaling system might also play roles in equilibrium reception (statocyst) or respiratory (the mantle and gills) functions in Aplysia 71 . It is expected that future studies will clarify these issues and provide a better understanding of the functions of the elevenin signaling system in protostomes. Finally, as elevenin-like peptides are present in cone snail venoms 66,67 , these studies might also provide insights that help promote the development of drugs.

Material and methods
Subjects and reagents. Aplysia californica (100-350 g) was obtained from Marinus, California, USA.
Aplysia are hermaphroditic (i.e., each animal has reproductive organs normally associated with both male and female sexes). Animals were maintained by circulating artificial seawater at 14-16 °C and the animal room was set with a 24 h light cycle with a light period from 6:00 am to 6:00 pm. All chemicals were purchased from Sigma-Aldrich unless otherwise stated.
Bioinformatic analysis of elevenin receptors. We first used NCBI to search for specific receptors and selected the nucleic acid sequences of receptors from species evolutionarily close to Aplysia, e.g., Platynereis. We performed BLASTn in NCBI and selected possible elevenin receptors according to E-value and similarity (positive). These sequences found by BLASTn were also searched in the AplysiaTools (http:// www. aplys iatoo ls. org/) to obtain more similar sequences for comparison. The complete nucleic acid sequences were searched for open reading frames (ORFs) in NCBI (https:// www. ncbi. nlm. nih. gov/ orffi nder/). Then the amino acid sequences of the putative receptors were subsequently analyzed in TMHMM-2.0 (https:// dtu. biolib. com/ DeepT MHMM) to determine whether they are a seven-transmembrane GPCR receptor. After the putative receptors were identified, phylogenetic trees were constructed using ClustalW alignment and the maximum likelihood method of 1000 replicates in MEGA X (https:// www. megas oftwa re. net/). Besides, we used the LG + G + F model to run all the phylogenetic trees. RNA extraction. The ganglia were removed from Aplysia, and Trizol was added to make up to 1 ml. The tissue was crushed by sonication, then, 200 μl chloroform was added. The solution was kept on ice for 15 min and was then centrifuged (10,000 × g, 4 °C, 15 min). The supernatant was transferred into a new centrifuge tube. An equal volume of pre-cooled isopropanol was added to the tube on ice. After standing at -20 °C overnight, the tube was centrifuged (10,000 × g, 4 °C, 10 min), and the supernatant was discarded. The pellet was washed twice with 1 ml of 75% ethanol/water and the centrifuge tube was shaken gently by hand to suspend the pellet. The tube was centrifuged again (12,000 × g, 4 °C, 10 min), the supernatant was discarded, and the precipitant was dried at room temperature for 5-10 min. Finally, 20 μl of nuclease-free water was added to dissolve the RNA and the concentration of RNA was measured.
Reverse transcription and PCR. RNA was reverse transcribed into cDNA using the reverse transcription kit (Takara, RR036A) according to the instructions. The cDNA obtained after reverse transcription was used as a template for PCR. Based on the nucleic acid sequence of the putative elevenin receptor, each pair of specific primers (see Supplementary Table 2) was designed in NCBI. We cloned the elevenin precursor and the three elevenin receptors. Both PCR reactions were performed at 98 °C/2 min pre-denaturation, 98 °C/10 s denaturation, ~ 60 °C (depending on specific primers: see Supplementary Table 2)/15 s annealing, 72 °C/30 s extension (depending on specific PCR products), and 72 °C/5 min for another 35 cycles. After running PCR, double digestion, ligation, plasmid transfection, plate plating, and single clone picking, the receptor was cloned into the vector pcDNA3.1( +). Then the elevenin receptor with pcDNA3.1( +) was sequenced and aligned to ensure the sequence is correct. www.nature.com/scientificreports/ IP1 accumulation assay. The metabolic inositol phosphate cascade is usually caused by the regulation of phospholipase C-β (PLC-β) associated with the Gαq subunit of heterotrimeric G proteins. Gαq-coupled GPCRs stimulate PLC-β activity, resulting in increased cellular D-Myo-Inositol 1-Phosphate (IP1). When the G proteincoupled receptor (GPCR) expressed in CHO-K1 cells (Procell, CL-0062) is activated, it is generated by the Gαq pathway through the appropriate ligand, and IP3 is formed which is degraded to IP1. Therefore, IP1 can be used as an alternative to IP3. First, the putative receptor cloned into the mammalian expression vector pcDNA3.1( +) of Aplysia was transiently expressed in CHO-K1 cells. CHO-K1 cells were cultured in F-12 K medium (Gibco, 21,127-022) containing 10% fetal bovine serum (Genial, G11-70,500) at 37 °C and 5% CO 2 . When the cell density in the 6-well plate reached 70-90%, 1.5 μg of putative receptor plasmid [in pcDNA3.1( +)] was mixed with 200 μl buffer (Polyplus) in a 1.5 ml EP tube and incubated at room temperature for 10 min. The DNA/ Transfection reagent (Polyplus) mixture was then added dropwise to the 6-well plate, mixed gently, and the cells were incubated overnight at 37 °C and 5% CO 2 . The next day, cells were digested, counted, diluted to a density of 40,0000 cells/ml, and replated in a white 384-well assay plate (Greiner bio-one, B14061RI), each well inoculated with the diluted liquid 50 μl and incubate overnight at 37 °C in 5% CO 2 . On the third day, the CHO cells were incubated for 1 h after adding peptide. The activation of putative receptors was detected by monitoring IP1 accumulation using the IP1 detection kit in Tecan Spark (Cisbio, 62IPAPEB). All procedures were performed according to the manufacturer's instructions for the IP1 assay kit. The elevenin peptide and analogs were synthesized by Sangon Biotech or Guoping Pharmaceutical (Supplementary Fig. 4) and were aliquoted in 50 nmol EP tubes and stored at -20 °C until use.
Statistical analysis. Dose-response curves and bar graphs for experimental data were plotted using Prism software (GraphPad). Data are expressed as the mean ± SEM. Final EC 50 values are rounded to two significant figures. All experimental data were taken from individual animals or preparations, and n refers to the number of preparations. Data were considered significant when P < 0.05. Statistical tests were performed using Prism software. They included Student's t-test (two groups) and one-way ANOVA (three or more groups). When oneway ANOVA shows a significant difference, individual comparisons with Bonferroni corrections were then performed.

Data availability
The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.